Mean field theory for double-perovskites 
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A mean field approximation of a model for double perovskites that takes into account the cou- 
pling between itinerant electron spins and localized spins is developed. As in previously reported 
theoretical results, and contrary to experimental observation, the critical temperature is suppressed 
for large electron density. An effective Heisenberg model reveals the cause of this discrepancy: the 
competition between degenerate antiferromagnetic and ferromagnetic channels. This degeneracy 
can be broken by the inclusion of a Hubbard-type U term. It is therefore suggested that electron 
correlation effects need to be incorporated in the minimal model of double perovskites in order to ex- 
plain the experimental observation of increasing ferromagnetic critical temperature with increasing 
electron doping. 

PACS numbers: 75.10.-b, 75.47.Gk 



I. INTRODUCTION 

A large magnetoresistance is a desirable property for 
certain technological applications such as writing and 
reading magnetic memories. Manganites, the most stud- 
ied colossal magnetoresistance systems^ have the draw- 
back of having low ferromagnetic transition (Curie) tem- 
peratures Tc (compared to room temperature). They 
are half-metals, namely, metallic for one spin orienta- 
tion and insulating for the other, but this property gets 
rapidly suppressed when increasing temperature. Half- 
metalicity produces the low field (extrinsic) magnetore- 
sistance measured in polycrystalline manganites and the 
tunneling magnetoresistance in artificially created bar- 
riers. Increasing the operation temperature of these 
devices is a major issue that has lead to the search 
of other half- metals with higher Tq. One such exam- 
ple are double perovskites of general formula A2BB'06 
(A=Sr,Ca,Ba,La,K, B=Fe, B'=Mo,Re)£2d£ Polycrys- 
talline double perovskites show large magnetoresistance 
at low fields due to half-metalicity even at room tem- 
perature, and their Tq is above 400K. 2 Band structure 
calculations2^*SiL£ reveal a ~ 0.8 eV gap in the major- 
ity up spin bands at the Fermi level while the down spin 
bands cross it A 

The ordered double perovskite structure consists of al- 
ternating BC>6 and B'Og octahedra in a cubic lattice. 
Anti-site disorder in this lattice has strong effects in mag- 
netic propertiesiSiiSiiiiiSiiii^ Fe is in the trivalent state 
(3d 5 ) with its five electrons localized in the spin- up d- 
orbitalsA Mo 5+ provides one electron (4c? 1 ) per Fe to 
the conduction band, electron density c = 1, and Re 5+ 
provides two (5d 2 ), c = 2. These electrons are mov- 
ing in the Fe-Mo/Re hybrid band formed by t2 g spin- 
down orbitals. The cubic symmetry causes a splitting 
between the t2 fl and e g orbitals, such that t2 S orbitals 
are lower in energy and the only ones occupied by the 
conduction electrons. The spin-down orbitals are well 
above the spin-up orbitals due to the strong Hund's cou- 



pling in Fe (see Fig. [I}. Fe is the magnetically active ion; 
its five localized electrons render a local spin S = 5/2. 
On the other hand, Mo and Re are paramagnetic ions^ 
The theoretical magnetization per formula unit is 
though smaller values have been measured probably due 
to anti-site disorder, as nearest neighbor Fe-Fe superex- 
change (SE) interaction is antiferromagnetic, 10,121 - 1 ^ In an 
ordered lattice, Fe-Fe SE is very weak as Fe ions are too 
far away (> 5.5A) for this interaction to be important. 
Doping on the A site changes the density of conduction 
electrons per Fe: substituting a trivalent ion (e.g. La) for 
Sr (S^-zLa^FeMoOg) increases the electron density to 
c = 1 + x while a monovalent ion (e.g. K) gives c = 1 — x. 

Superexchange interaction, which is extremely weak 
in double perovskites, has been ruled out as the cause 
of magnetic ordering in these materials. Instead, there 
is experimental evidence for the existence of two sub- 
lattices^ the localized spins in Fe and the delocalized 
electrons, that interact antiferromagnetically due to the 
strong Hund's coupling J on Fe ions. Other materi- 
als that present these two coupled sublattices are man- 
ganitesi and diluted magnetic semiconductors such as 
Gai^^Mn^AsiiSiii In manganites the conduction elec- 
trons are ferromagnetically coupled to the localized Mn 
spins that form a pseudocubic lattice. Hund's coupling 
in manganites is very large and, in the limit J — ► oo, 
leads to the double exchange (DE) mechanism produc- 
ing ferromagnetism. DE basically implies that the spin 
of the conduction electron follows the orientation of the 
local spin. As hopping does not flip the spin, the kinetic 
energy of the conduction electrons is minimized when the 
local spins are all parallel to each other* 18 Hund's cou- 
pling in double perovskites is antiferromagnetic^ as the 
spin-up Fe levels are fully occupied, but it can equally 
lead to ferromagnetic order. However, unlike in man- 
ganites, the strong Hund's coupling only applies to every 
other site, since Mo is paramagnetic^ and, therefore, the 
minimal model for double-perovskites is different from 
the simple DE in manganites^ 



Susceptibility measurements of double-perovskites in 
the paramagnetic regime give a positive Curie- Weiss 
parameter (i.e. ferromagnetic interaction) piLiSiiS This 
rules out a superexchange scenario, since SE would only 
give ferromagnetic Fe ordering through antiferromagnetic 
coupling between Fe:3d 5 and MoiSd 1 (ferrimagnetic or- 
dering) and in this case a negative Curie- Weiss parame- 
ter should be observed. A conduction electron mediated 
ferromagnetism is suggested by the observation that the 
strength of the ferromagnetic coupling and thus, the To, 
can be increased by electron dopingii&SfliSAiS In turn, the 
increase in Tq has been shown to be accompanied by an 
increase of the density of states at the Fermi levels 

A model for double perovskites that takes into ac- 
count the coupling between itinerant electrons and local- 
ized spins on Fe has previously been studied theoretically 
by means of dynamical mean-field theory^ and Monte- 
Carlo simulations^ Their common result is that the crit- 
ical temperature Tq is suppressed as the electron density 
increases above a certain value. This is clearly in contra- 
diction with experimental observations of increasing Tq 
by electron doping. On the other hand, similar models 
for perovskite manganites and diluted magnetic semicon- 
ductors (GaMnAs) find that the Tq increases with the 
density of states. 
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FIG. 1: (color online) Description of the parameters of the 
model and definition of the paramagnetic state used, t stands 
for tF C -Mo- The level with energy — J — A correspond to the 
localized electrons in _Fe(3d 5 ) that are considered in the model 
as localized classical spins S = 5/2. The dash box defines the 
unit cell i. In the limit J — ► oo the spin of the conduction 
electron on the upper band is always strictly antiparallel to 
the localized spin. 



In this paper a mean field theory approximation is de- 
veloped on this model for ordered double perovskites. 
The previously published results^^ are recovered in the 
appropriate limits. As mentioned already, these theoret- 
ical results are not consistent with experimental data. 
The advantage of the mean field theory developed here 
is that the reason for the failure of this model becomes 
transparent. The Hamiltonian can be written as an ef- 
fective Heisenberg model with ferromagnetic and anti- 
ferromagnetic terms. Due to the degeneracy of the spin 
up and down levels in the paramagnetic atom, these two 
channels compete resulting in the suppression of Tq for 
large enough values of the electron density. We show 
how the intraband Coulomb repulsion U, which penal- 
izes the occupation of two spin orientations at the same 
site, can prefer the FM channel over the AF leading to an 
increase of Tq with doping, as observed experimentally. 
This work will not take into account the possibility that 
the Coulomb repulsion also induces orbital order m& 

This paper is organized as follows. In Sec. II the min- 
imal model for double perovskites is introduced. In Sec. 
Ill, the mean field approximation is described, the model 
is written as an effective Heisenberg model, and the ef- 
fects of adding electron correlation effects are analyzed. 
In Sec. IV, an alternative mean field calculation is de- 
scribed. This approach integrates out the effect of the 
paramagnetic sites and considers only the Fe sites, mak- 
ing a clear connection with the double exchange model 
used in the description of manganites. We conclude in 
Sec. V. 



II. MODEL 

The full Hamiltonian can be written as the sum of two 
terms 

H = i?KE + ^on-site- (1) 

The first term is the kinetic energy of the conduction 
electrons. The second term takes account of the strong 
Hund's coupling on Fe sites, J, and the difference in elec- 
tronegativities, A, between Fe and Mo/Re sites. 

The conduction electrons move between t2 S orbitals. 
These orbitals have planar symmetry which implies hy- 
bridization only takes place between t2 S orbitals of the 
same symmetry. Therefore, the kinetic energy consists 
of three degenerate, two dimensional, tight binding sys- 
tems xy, yz, and zx, H KE = Hf® + Hf® + Hf^. The 
matrix element iF C -Mo connects the d a b (a, b = x,y,z) 
orbitals of nearest neighbors in the ab plane. irvio-Mo 
connects nearest neighbors in the Mo sublattice. t^e-Fe 
is expected to be very small due to the more localized 
nature of the 3d states, and is neglected here. 

The large local spins on Fe sites (S = 5/2) can be 
considered classical and are characterized by an angle 9i. 
Due to the large value for the Hund's coupling, the spin 
of the conduction electron on Fe follows adiabatically the 
classical local spin configuration. 
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III. MEAN FIELD APPROXIMATION. 
DISPERSIVE MO BANDS. 

Tq is calculated from the expansion of the Free Energy 
F in powers of the magnetization m = (cos0j), which is 
very small close to Tq. Tq is defined by the condition 
d 2 F/dm 2 = 0. In order to write the Free Energy, we 
have to define the paramagnetic regime correctly. Diago- 
nalization of the one-dimensional problem shows that the 
density of states has a gap in the paramagnetic regime. 
This gap is related to the fact that Mo is paramagnetic 
and there are two conducting channels: one through each 
spin orientation. 

In the paramagnetic regime, the spins on the Fe be- 
come randomly oriented. Mo is paramagnetic, allowing 
two spin states. To recover the gap in the DOS, the orien- 
tation of these two states in the Mo has to be referenced 
to a neighboring Fe spin, as shown in Fig.^ The Fe and 



Mo related in this way constitute our unit cell. In this 
way, the spin of the electrons in the Mo is chosen to be 
parallel or antiparallel to the direction of the Fe spin in 
the same cell. Inside the unit cell, the Fe-Mo hopping is 
1 (for parallel spins) or (for anti-parallel spins) and be- 
tween different cells it is determined by the angle formed 
by Fe ions in neighboring cells 6^ = 9i — 6f cos(9ij/2) 
(spin t channel) and sin(0jj/2) (spin j channel). There- 
fore, the system consists of two different channels that 
could not be distinguished if the relative orientation of 
the spins in Mo and Fe were not correctly defined. In 
the virtual crystal approximation, each site sees an av- 
erage of all the sites in the lattice so the relevant coef- 
ficients of the hopping terms are the thermal averages 
(cos(%/2)) = (cos(0/2)) and (sin(%/2)) = (sin(0/2)). 
In the paramagnetic regime both are equal to 2/3. 

For each of the equivalent planes, the Hamiltonian can 
then be written 



H xy = ( J - A) ^ d t d i + 'Mo-Mo X] \ C0S 4) ( C i,p C i'P + C tap C j,o.p) ~ ( Sin ^/ ( C i,ap C 3,P + C t,p C j,ip) 



^Fo-Mo Y^ d i 



•:.3> 

^i,p ~t~ ( COS — y {Ci— x ,p H~ Ci — x — y,p Ci — y,p) ( Sill — y \^i—x,ap ^i—x—y.ap H~ Ci— y,ap) 



,(2) 



where c% v ( a p\ destroys an electron in Mo at cell i with 
spin parallel (antiparallel) to the spin of the Fe core spin, 
and di destroys an electron in Fe at cell i, with the spin 
parallel to the core spin. H yz and H zx have identical 
form and give the same contribution to the total energy. 

Close to the magnetic transition, we can write the hop- 
ping coefficients (cos(0/2)) and (sin(0/2)) as an expan- 
sion in m to second order 




(3) 



where we are following the same procedure as in 
Ref. Using these expressions in the Hamiltonian 

in Eq. (J2J, and taking into account that there are three 
equivalent 2-dimensional bands, the kinetic energy can 
be calculated as a function of m: Eke = -E'ke + X m2 ■ 
Knowing that the entropy of the spin system is 

where h is an external magnetic field and /3 = k^Tc, and 
that d 2 F/dm 2 = 0, we get T c = 2/3 [dE KE /d{m 2 )} [H 
|29| (see Appendix . 

Numerical results for Tc with parameters J — A = 0.3 
eV, £mo-Mo = 0.15 eV, and iF C _Mo — 0.3 eV are shown 



in Fig. |21 These parameters are consistent with ab-initio 
calculations^ and similar to the ones used in previous 
theoretical works on this modeli^S J + A is consid- 
ered to be infinite as the transitions to the Fe spin level 
parallel to the localized spin involve very large energies^ 
The results are in agreement with Monte-CarloSt and dy- 
namical mean-field theory^ calculations but in disagree- 
ment with experiments. In general, the calculated Tc 
is lower than that measured in experiments and shows 
a maximum around c *~ 1, while it is suppressed for 
larger electron densities c ~ 2. This behavior persists 
for a wide range of parameters (0 < <Mo-Mo < 0.25 eV, 
< i Fc -Mo < 0.5 eV, and -1.5eV < J - A < 1.5eV). 

The maximum of Tc around c = 1 and its suppression 
around c = 2 is related to the form of the density of 
states. Around c = 1 there is a maximum on the DOS 
and, close to c = 2 both parallel and anti-parallel spin 
bands from Mo are filled. The effect of this filling is more 
easily understood by introducing an effective Heisenberg 
model. 



A. Effective Heisenberg model 

To analyze the source of the discrepancy between the- 
ory and experiment the energy is written as a function 
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FIG. 2: (color online) Tc versus electron density c for a par- 
ticular set of parameters J — A = 0.3eV, tMo-Mo = 0.15eV, 
and tF C _Mo = 0.3eV. For all the range of parameters studied, 
there is a peak in Tc around c ~ 1 and Tc = at c ~ 2. 
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where the J's are the expectation values of the operator 
pairs in Eq. J5J. In the m — > limit, Eq. JSJ is an effective 
Heisenberg model, 
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(6) 



Therefore, this effective Heisenberg model has com- 
peting ferromagnetic (Jc's) and antiferromagnetic (Js's) 
terms. In Fig. Eta) the values of the four different cou- 
plings are plotted as a function of electron density. As 
shown in Fig. GJb), the total coupling J e g is antiferro- 
magnetic for large values of electron density. The elec- 
tron density at which J e g becomes zero does not change 
significantly within the range of the tight-binding param- 
eters used. Therefore, we cannot expect to override the 
suppression of Tq within this model. The competition 
between the ferromagnetic and antiferromagnetic chan- 
nels can lead to phase separation^ due to the fact that 
J e ff (and, consequently, Tc) depend strongly on the elec- 
tron density^ 

The ferromagnetic and antiferromagnetic channels are 
degenerate in energy since Mo is paramagnetic and both 
parallel and antiparallel states are equally populated. 



FIG. 3: (color online) (a) Heisenberg coupling values for the 
ferromagnetic (solid) and antiferromagnetic (dashed) chan- 
nels. The antiferromagnetic coupling gets stronger for large 
electron density c > 1.5. (b) The total coupling J e ff as defined 
in Eq. Q. 



be broken to obtain agreement with experimental results. 
Therefore, double occupancy of a site must be penalized 
which is, of course, expected in the presence of Coulomb 
interaction induced on-site electron correlations (so far 
neglected in the theory). This can be done with an in- 
traband Hubbard term of the form 



(8) 



Us- 



In Fig. 21 is plotted for different values of U. 
ing relatively small values for U (U << W ~ 8t), we 
observe that Tc is strongly enhanced and the position of 
the maximum is shifted up from c ~ 1, consistent with 
experiment. 



IV. MEAN FIELD APPROXIMATION. NON 
DISPERSIVE MO BANDS. 

If we neglect the direct hopping between the Mo or Re 
orbitals, the connection between the model studied here 
and the double exchange model used in the description 
of the manganites becomes transparent. The Mo, Re 
orbitals can be replaced by an energy dependent direct 
Fe-Fe hopping, and a correction to the energy of the Fe 
orbitals. These quantities are: 



B. Coulomb interaction 

It is clear from the previous section that the degeneracy 
of ferromagnetic and antiferromagnetic channels needs to 
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FIG. 4: (color online) Tc versus electron density c for different 
values of the Hubbard U. J — A = 0.3 eV, tMo-Mo = 0.15 eV, 
and tF C _Mo = 0.3 eV. Moderate values of Coulomb repulsion 
U produce a large increase in critical temperature and its 
suppression is shifted to higher electron density. 



When the separation between the Fe and Mo levels, J — 
A, is large compared to the effective hybridization of the 
Fe levels, namely, 



IJ-AI « 



''Fc-Mo 
J — A\ ' 



(10) 



the model reduces to an effective double exchange 
model, 27 with hopping integral: 



*Fe-Mo COs(%/2) 

U-AI 



(11) 



A mean field solution to the problem can be obtained by 
making the substitution: 



cos 



(12) 



where m is the temperature dependent magnetization of 
the Fe ions. 

Keeping the full energy dependence of the tight bind- 
ing parameters in Eq. © , one finds a self consistent so- 
lution for the band structure, which can be written as 
a quadratic equation for the band energies, ■ . This 
equation can be solved, and we obtain: 



A- J 

k X ky ~ J + 2 i 



(A-J)= 



4f 2 

^ Fc-Mo 



1 + 



-1 + (cos(fe x ) + cos(k y )Y 



(13) 



These equations give the exact solution of the model at 
zero temperature (m = 1), and they should describe qual- 
itatively the changes in the electronic structure induced 
by the magnetic fluctuations. Note that the absence of 
direct hopping between Mo orbitals leads to a dispersion- 
less band which has only weight at the Mo sites. 



From the knowledge of the dependence of the electronic 
bands as function of the magnetization, one can calculate 
the electronic contribution to the free energy, and obtain 
the value of the Curie temperature, as discussed in Ap- 
pendix A. We find: 



Tc 



2 d 2 E KE 

3 d(m) 2 



I /< **£ 

J Je k <e F ■ 



="Fc-M 



|[cos(A; x ) + cos(fc y )] 2 - 1 j 



(A-./) 2 



z(, Fc-Mo 



o {l + \[l [-1 + (cos(fc x ) + cos(fc y )) 2 ] } 



(14) 



The main drawback of this approximation is that it 
does not allow to estimate the contribution of the direct 
hopping between Mo orbitals, iMo-Mo, which is expected 
to be comparable to ipe-Mo- These terms tend to reduce 



the dependence of the electronic energy on the magne- 
tization of the Fe sites, lowering the value of Tc- The 
qualitative dependence of Tc on band filling obtained in 
this approximation is, however, similar to the one con- 
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sidcred in Sec. IIIII with Tq — ► for electronic densities 
c rs 2. This fact can be explained, qualitatively, by not- 
ing that the value of Tc in the effective double exchange 
model derived here tends to have a maximum when the 
band is half filled. Because of the triple degeneracy of 
the Mo orbitals, the corresponding density is c = 1.5. 



V. CONCLUSION 

We have developed a mean field theory for double per- 
ovskites, e.g. S^-^LazFeMoOe, within a minimal ef- 
fective model (Fig. ^| including the strong Hund's cou- 
pling on the Fe sites and the various contributions to 
the kinetic energy of electron hopping between t2 g or- 
bitals through the Mo sites. Ferromagnetism arises in 
the system due to constraints imposed on the hopping 
kinetic energy, rather than due to Fe-Fe superexchange. 
Our simple mean field theory reproduces the earlier the- 
oretical results obtained with dynamical mean field ap- 
proximation^ and direct Monte Carlo simulations. 25 Our 
theory gives a reasonable semi-quantitative description of 
the observed Curie temperature in the double perovskites 
for the 'undoped' x — system where the carrier density 
c = 1 (per unit cell). But, for the doped double per- 
ovskites (x ^ 0) our mean field theory, along with the 
existing theories of Refs. 1 2 'll and l25l predicts a decreasing 
Tq with increasing x (with c = 1 + x for electron dop- 
ing and c = 1 — x for hole doping) , which disagrees with 
experimental observations. The experimental finding is 
that Tq increases with doping whereas the theory finds 
a maximum around c ~ 1 (i.e. x = 0). 

We suggest, based in our mean-field formalism, that 
the experimental observation of increasing Tc with dop- 
ing may be an electron correlation effect which opposes 
double occupancy of sites due to intraband Coulomb re- 
pulsion. By introducing a simple Hubbard-U type in- 
traband correlation term, we qualitatively reproduce the 
experimental trend of increasing Curie temperature with 
increasing doping. In addition, the introduction of the 
Hubbard-U term also enhances the Tq itself bringing the- 
ory and experiment into better quantitative agreement. 
We therefore believe that strong correlation effects are 
an inherent property of double perovskites. 

Finally, we discuss the possible connection between 
double perovskites (DP) and diluted magnetic semicon- 
ductors (DMS), e.g. Gai^Mn^As, from the perspec- 
tive of magnetism. At first, one notices some superficial 
similarities between DP and DMS materials: both have 
optimal Tq of the order of a few hundreds of Kelvin (al- 
though the highest reported Tq in GaMnAs is around 
200K, substantially below the room temperature 2c 's 
routinely seen in DP materials); both manifest Tc's in- 
creasing with doping, thereby indicating a role for carrier 
mediated ferromagnetism. There are, however, impor- 
tant differences between DMS and DP magnetic proper- 
ties. In DMS, Tq — for x = since Mn atoms serve 
the dual roles of dopants (providing the carriers, which 



are holes for Gai-^Mn^As) and magnetic moments (i.e. 
the long range ferromagnetic order arises from the order 
of the local Mn moments), and therefore ferromagnetism 
vanishes in the absence of Mn. Thus the x — situ- 
ation in DMS is qualitatively similar to the x = 1 DP 
situation. The common modeUSiii adopted in the liter- 
ature to understand DMS ferromagnetism is a carrier- 
mediated RKKY-Zener indirect exchange coupling be- 
tween the Mn moments, with the mean-field DMS Tq 
given by Tc ~ | Jpd^xn 1 / 3 , where J p a is the local 'pd' 
exchange coupling between the Mn d-level and the p- 
type holes (with density n) in the valence band of GaAs. 
This RKKY-Zener type DMS mean-field theory is obvi- 
ously completely inapplicable to the DP materials as it 
would predict an absurd DP Tc of 10 5 — 10 6 K or larger 
(since both the magnetic moment density and the car- 
rier density are substantially higher in DP materials than 
in DMS materials). We have developed the appropriate 
DP mean-field theory in this paper with a reasonable 
T c - 10 2 - 10 3 K. 

It has, in fact, been suggested in the literature 3 ^ that 
the DMS systems may actually be closer to the non- 
perturbative double-exchange limit that the perturba- 
tive RKKY limit. In such case, the DMS and the DP 
systems are more similar in nature (and they are both 
then closer to manganites, which are the quintessential 
double-exchange materials^i* 3 ^) . But, in this limit, in- 
creasing doping should invariably lead to the eventual 
suppression of Tc, as we find in the theory developed 
in this paper. An important difference between DP and 
DMS materials is, however, the fact that the DMS sys- 
tems lose their ferromagnetism (i.e. Tc becomes zero) 
for large values (x ~ 0.1) of Mn concentration. This 
also sharply distinguishes the DMS and the DP materi- 
als. A natural question, based on our argument in favor 
of the important role of a Hubbard type U-term in the 
DP materials is whether such electron correlation effects 
are also important in DMS materials. The answer to 
this question is not obvious at this stage. One possi- 
bility is that correlation effects are completely negligible 
in the DMS systems since the effective carrier density 
is extremely low (n ~ 10 19 — 10 20 cm~ 3 in DMS com- 
pared with 10 22 — 10 23 cm~ 3 in DP materials), making 
the physics of double occupancy irrelevant. Much more 
work will obviously be needed to further understand the 
relationships and the differences in the magnetic prop- 
erties and mechanisms for various 'oxide-type' magnetic 
materials such as manganites, double perovskites, diluted 
magnetic semiconductors, magnetically doped magnetic 
oxides (e.g. Tii_ a; Co 2 ;02), and even systems such as 
Fei-zCo^Si where correlation effects are thought to play 
an important role. 
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APPENDIX A: MEAN FIELD EXPRESSION FOR 
THE CURIE TEMPERATURE 



where %m 2 is the kinetic energy of the conduction elec- 
trons. 

To obtain h, we minimize the total free energy with 
respect to h, dF total /dh = 0. In the limit of small h 



The Free Energy of a system of classical spins of mag- 
nitude unity in an external magnetic field h is 



1 / sint 
F= -/? l0S ( 2 -/^ J' 



(Al) 



from where the magnetization can be calculated as 
dF 1 1 

m =< m >= 



dh tanh(/3/i) 0h 
The entropy of the spin system is then 



(A2) 



mj3h 



(A3) 

The total energy of the system, assuming the electrons 
are at zero temperature and neglecting direct interactions 
between Fe spins, can be written 



F 



total 



0h J 



mfth 
(A4) 



sinh(/?/i) = (/3h) 2 
0h = + 6 



(A5) 



In this limit, the minimization condition gives f3h 
—3m, and the free energy gets the form 



F 



total 



-Ilog(2) + m 2 [x+\k B T 



(A6) 



Therefore, the Curie temperature i 



,29 



(A7) 
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